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Abstract 

An improvement in the treatment of the isovector channel of relativistic mean field (RMF) 
models based on effective field theory (E-RMF) is suggested, by adding an isovector scalar ( 6 ) 
meson and using a similar procedure to the one used by Horowitz and Piekarewicz to adjust the 
isovector-vector channel in order to achieve a softer density dependent symmetry energy of the 
nuclear matter at high density. Their effects on the equation of state (EOS) at high density and 
on the neutrino mean free path (NMFP) in neutron stars are discussed. 
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I. INTRODUCTION 


Neutrino transport in stellar matter plays an important role in some phenomena such 
as the mechanism of supernovae explosions, structure of protoneutron stars, etc. Theoret¬ 
ical input needed for understanding the neutrino transport is the neutrino mean free path 
(NMFP) and the equation of state (EOS). A unified matter model used in calculating both 
observables is a requirement for having a consistent neutrino transport result. We note that 
considerable efforts have been devoted to investigate the neutrino interaction in matter at 
high density 2, 3, 4, 3, 6, 7,1a, il llli III 12, 13, Im llf| with different matter models, variety 
of approximation levels and purposes. 

To this end, a kind of matter models that can be used is the relativistic mean held 
(RMF) model. The original or standard RMF models use the sigma, omega, and rho mesons 
with additional cubic and quartic nonlinearities of sigma meson to effectively describe the 


interaction among nucleons. NL-3 parameter set m belongs to this model. The parameter 
set has been very successful in the description of a variety of ground state properties of 
nuclei [Q. Details of the models as well as their applications can be seen in Refs. IihI. liil I 20 I. 
The effects of including delta meson in standard RMF models including the corresponding 
linear responses have been studied for the asymmetric nuclear matter at low densities in 


Refs. |2!|, !22j, |23j, whereas for heavy ion collisions in Refs. 


M 


25 


26] and on the neutron 


star properties in Ref. [27]. It was found that the 5 held leads to a larger repulsion in the 
dense neutron rich matter (stiffer symmetry energy that leads to a larger proton fraction at 
high density), as well as a definite splitting of proton and neutron effective masses. Both 


features are influencing the stability conditions of a neutron star 


27]- 


Inspired by the concept of effective held and density functional theories for hadrons, 


Furnstahl, Serot and Tang 


28| constructed a new RMF model (from now on, will be denoted 


by E-RMF) that can be considered as an extension of the standard RMF models. One of 
the parameter sets in this model is G2. Besides yielding accurate predictions in finite nuclei 


and normal nuclear matter 
of quartic sigma meson coup 


12! QQ, 


G2 has the interesting features like a positive value 


spectrum of this model 


ing constant that leads to the existing lower bound in energy 


30 3: 

t 3 


31] and to the missing zero sound mode in the high density 
Moreover, the agreement of the nuclear matter and the 


symmetric nuclear matter 
neutron matter EOS at high density of G2 with the Dirac Brueckner Hartree Fock (DBHF) 
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calculation 1211I30l] is better than those of NL-3, NL-1 and TM1 models (the standard RMF 
plus a quartic omega meson interaction). Nevertheless, from the comparison between the 
neutron matter EOS of this model and that of the DBHF result, the authors of Ref. 2p| 
pointed out that the present form of the E-RMF model still needs a substantial improvement 
in the treatment of the isovector sector. It has also been shown that the G2 parameters 


set of the E-RMF model still predicts a too large proton frac 


ion 


B. 


It is known that 


proton fraction correlates to the direct URCA cooling process 133] . It is also known that 
this problem is caused by the role of isovector terms. So far, the effects of the delta meson 
inclusion on this model has not been studied yet. 

Therefore, in this work, first, we will study the effects of adding an isovector-scalar (5) 
meson in the E-RMF model and the effects of the isovector-vector channel adjustment by 


using a similar procedure to the one used in Ref. [1]. The aim of these adjustments is to 
achieve a softer density dependent symmetry energy of nuclear matter at high density. The 
symmetry energy has a wide range of effects, such as from giant dipole resonances to heavy 
ion collisions in nuclear physics and from supernovae to neutron stars in astrophysics. In spite 
of its diverse influences, its magnitude and density dependence are not well understood |2J]. 
More detail informations on the role of symmetry energy and related topics can be consulted 
to Ref. 34] and references therein. Second, we will also extend the analysis of our previous 
report [4] to give a more solid argument about the source of the possible appearance of the 
anomalous behavior in the NMFP of neutron stars predicted by RMF models. Here the 
anomalous behavior in the NMFP means a contra intuitive NMFP results in a form of the 


0 


decreasing of the mean free path with respect to the decreasing of the matter density 
It has been known that this anomalous behavior exists in the NMFP predicted by non 
relativistic nuclear models of the Skyrme type j^j, [j| . This anomalous behavior is attributed 
to the dominating term at high density that responsible to the appearance of the acausal 
behavior (the speed of sound exceeds the speed of light) of the model at high densities 


It has been reported in Refs. 


0 


that relativistic models alleviate this problem. But we 


have eventually found that not every parameter set of relativistic models is free from this 


Q 


problem 

In section |H| we will briefly explain the formalism used in this work. In section IIIII we 
discuss the results of our calculations. We will give the summary of our findings in section llVl 
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II. FORMALISM 


This section contains a very brief description of the self consistent models for nucleons 
plus the standard non-interacting Lagrangian densities for electron and muon as well as the 


;he standard model of weak 
is used, i.e., the ground 


interaction between neutrino-electron with matter based on 
interaction. Here, a similar assumption with that of Refs. 3, 
state of the neutron star is reached once the temperature has decreased below a few MeV. 
This state is gradually reached from the later stages of the cooling phase. The system is 
then quite dense and cool so that the zero temperature approximation is valid. In this 
case the direct URCA neutrino-neutron scattering is kinematically possible for low energy 
neutrinos at and above the threshold density when the proton fraction exceeds 1/9 Q or 
slightly larger if muons are present. Furthermore, the absorption reaction is suppressed. For 
simplicity, we neglect the RPA correlations. 


A. E-RMF Model 

The calculation is done in the framework of the relativistic mean held approximation, 
effective Lagrangian density used to describe nucleons interactions is taken from Refs. 


I2S, 


35]. This Lagrangian is constructed with a nonlinear realization of chiral symmetry J28(. 


The explicit form of the E-RMF effective Lagrangian density reads 


£ nuc = C N + C M . 


( 1 ) 


For nucleons, the Lagrangian density up to order v = 3 is given by 

C N = + ivp + ig p b p + igJC^) + gja^^a^ - M 

+ 9a(r}^ - 


where 


=(*’:) > v n = 


( 2 ) 


(3) 




£ = exp(in(x)/fn), tt(x) = -f • n(x), 


(4) 

(5) 
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( 6 ) 


k(x) = f' A x )i 


b P u D p b u T ig p [b p , b j,], .D^ <9^ -I- 777 ^, 


(7) 


^ = d p V u 


d„V„, 


( 8 ) 


CT 


/il' _ 


“sir, 71 - 


(9) 


Here, p, n and M are the proton-, neutron-field and nucleon mass, while cr, if, and 
are the cr, 7 r, u; and p meson fields, respectively. For mesons, the Lagrangian density up to 
order v = 4 reads 

£„ = iyTV^t/cW*) + i/pHf/C/t - 2) + iS^^«r - ^(5^) - jV (u V<“' 


- TV(5^^) + 4(1 + + f f)*l" + jCosifVC") 2 


m 


+ P + - ™> 2 a +f ^ + =f^>, 


( 10 ) 


where 


77 c^zy, du^fi T ^ly[xi Vp\ 7[o^, o^]. (11) 

In the mean field approximation, the it meson does not have a contribution. If we set 771 , 


772 , Coi Vp and /p equal to zero, in the Lagrangian density, the same equations of state for 
nucleons and mesons of the standard RMF models [la, [ljj, |20[ can be obtained. 

The density dependence of the modified nuclear matter symmetry energy after including 
the isovector-vector nonlinear term, 


£hp = 4A vglgl w ■ K R'%, 


7 pjco 


( 12 ) 


in the Lagrangian density of the standard RMF model and then followed by an adjustment of 
g p and Ay parameters has been for the first time calculated by Horowitz and Piekarewicz |l|. 
Motivated by a similar philosophy, other additional nonlinear terms, but with different forms, 


have also 
in Refs. |1 


oeen studied in Ref. [2f|. In this paper, we follow the same procedure as given 
, 13(1, but since we use the E-RMF model which already contains an isovector- 


vector nonlinear term, the density dependence of the nuclear matter symmetry energy can be 
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modified without adding a new nonlinear term, instead, it only needs an adjustment of the 
g p and rj p parameters while keeping the symmetry energy at the same value, i.e., E sym = 24.1 
MeV at kp =1.14 fm -1 . The argument behind this procedure is that the symmetry energy 
at the saturation density (kp =1.32 fm -1 ) is not well constrained experimentally. However, 
an average of symmetry energy at full density (the average density is less than saturation 


density) and at surface symmetry energy is constrained by binding energy of nuclei [1|, |3 
To study the effects of a <5 meson addition in the E-RMF model, we add to the Lagrangian 
density of that model the following terms, 


£s = (?) +'tpg S T-S'i/j. 


For electron and muon, the Lagrangian density reads 

Y -rni)l. 

l=e~, 


(13) 


(14) 


All matter properties used in this work can be derived from these Lagrangian densities. A 


more detail step of the derivations can be seen in Refs. 


1 28, Q 30, 35, 36] 


Coupling 


constants for all parameter sets used in this work are displayed in Table 0 To determine the 
fraction of every constituent, we imposed the requirement that in the neutron star at zero 
temperature, the chemical potential is in equilibrium and the charge is neutral. 


B. Neutrino Mean Free Path 

For neutrino-electron matter interactions based on the standard model of weak interac¬ 
tion, the Lagrangian density for every constituent is 

£int = ^=[ p 7 M (! - ( 15 ) 

where J ^ = 7 p (Cy — C\ 75 ) and j— n,p, e~,n~. The values of Cy and C\ can be seen in 
Table EU 

The neutrino differential scattering cross-section can be derived from the Lagrangian 
density and has the form of 


1 d 3 a 
V d?Vt'dE' u 


G f E' u 
32vr 2 E v 


im(vn^). 


(16) 
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TABLE I: Coupling constants of the parameter sets used in this work. 


Parameter 

G2 

NL-3 

NL-Z 

Z271 

ms/M 

0.554 

0.541 

0.520 

0.495 

9s/( 4tt) 

0.835 

0.813 

0.801 

0.560 

flv/(4 tt) 

1.016 

1.024 

1.028 

0.670 

9r/ ( 4? r) 

0.755 

0.712 

0.771 

0.792 

« 3 

3.247 

1.465 

2.084 

1.325 

Av4 

0.632 

-5.668 

-8.804 

31.522 

Co 

2.642 

0 

0 

4.241 

m 

0.650 

0 

0 

0 

m 

0.110 

0 

0 

0 

v P 

0.390 

0 

0 

0 


TABLE II: Coupling 
0.223 and gA = 1.260 


constants of neutrino-electron matter interactions. 



Here we use sin 2 9 W = 


Target 

Cy 

C A 

n 

-0.5 

-9a/2 

P 

0.5 - 2 sin 2 6 W 

9a/2 

e 

0.5 + 2 sin 2 9 W 

1/2 

P 

—0.5 + 2 sin 2 9 W 

-1/2 


Here E u and E' v are the initial and final neutrino-electron energies, respectively, Gp = 
1.023x10 _5 /M 2 is the weak coupling, and M is the nucleon mass. The neutrino-electron 
tensor can be written as 

L^v = 8 [ 2 kftku + {k.q)g,j, v - ( k + q^k v ) 

~ ( 17 ) 

where k is the initial neutrino-electron four-momentum and q — (qo, q) is the four-momentum 
transfer. The polarization tensor n M ", which defines the target particles, can be written as 

n U(q) = -if -SEi V[G>'(p)J'G>'(p + 9 ) 4 ], (18) 
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where p = (po,p) is the corresponding initial four-momentum, and G(p) is the target particle 
propagator. The explicit form of the nucleon propagator is 

1 


G n ’ p (p ) = {f + M* 


l 7T 

- + — 


p* 2 - M* 2 + ie E 


x 6(p‘„-E*) 0 (j% n - |p|) 


(19) 


where E* = E + So is the nucleon effective energy and M* = E + Eg is the nucleon effective 
mass. The S 0 and Eg are the scalar and time like self energies, respectively. Electron 
and muon propagators have similar expressions, except the effective (starred) quantities in 
Eq. (fTTfl) are replaced by the free ones. 

The NMFP (symbolized by A) as a function of the initial neutrino energy at a certain 
density is obtained by integrating the cross section over the time- and vector-component of 
the neutrino momentum transfer, as [t| 

rZEu-qo r 2 E„ |^| -\ 


A (E v 


d\q\ 


2Ev \q\ 1 

dq 0 ^Mz-2TT 


'90 


'0 E' V E„ 


V d 2 9JdE' 


( 20 ) 


Equation EH) is used to calculate the NMFP in the neutron star matter of Figs. 0E1 
andd 

III. RESULTS AND DISCUSSIONS 


In this section we study the effects of the isovector-vector channel adjustment and the 
addition of 6 meson in E-RMF models on the corresponding nuclear matter properties pre¬ 
dictions. We also study the role of every factor involved in the predicted NMFP in neutron 
stars of RMF models. We start with considering the effects of isovector-vector channel of 
the E-RMF model. Here we analyze the effects of different g p and rj p combinations of the G2 
parameter set. The value of various coupling constant combinations can be seen in Table Hill 
Their effects on matter properties are shown in Fig. |T1 

The symmetry energies of sets I-IV given in Table UTTI are shown in the upper left panel. 
Since it has been pointed out by Lattimer et al. |23], that the crucial role of the proton 
fraction value for the onset of the direct URCA process is to enhance the neutron star 

n 

cooling rate {42|, we plot the corresponding proton fractions in the upper right panel. To 
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TABLE III: Isovector-vector channel adjustment in the G2 parameter set of E-RMF model. 


Isovector 



Set 


parameter 

I 

II 

III 

IV 

9 P 

9.358 

9.483 

11.786 

13.687 

hp 

0.190 

0.390 

4.490 

8.490 




/Wai 


/Wa) 



e(GeV fm 3 ) 



Pb^Po 


FIG. 1: Isovector-vector channel adjustment in the E-RMF model. Effects of different combina¬ 
tions of g p and r] p on the symmetry energy of the nuclear matter are shown in the upper left panel, 
on the pressure and M* of the PNM are in the lower left and right panels, respectively, and on 
the proton fraction predictions in the upper right panel. Shaded region in the upper right panel 
corresponds to the proton fraction threshold of the direct URCA process Q. 

estimate their effects on the EOS of neutron star, we plot the pressure as a function of 
energy of the pure neutron matter (PNM) as the dominant contributor in the EOS of a 
neutron star, in the lower left panel. It is clear from the figure in the lower right panel that 
the adjustments of g p and g p have no effect on the PNM effective mass (M*). 
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TABLE IV: Isovector-vector channel adjustment in the Z271 parameter set of the Horowitz and 
Piekarewicz model ll. 


Isovector 


Set 


parameter j 

II 

III 

IV 

9 P 9.498 

9.672 

11.506 

12.145 

Ay 0 

0.01 

0.03 

0.035 


In conclusion, for the E-RMF model, proton fraction predictions depend strongly on the 
behavior of the density dependent of E s ym at high density. On the other hand, the predicted 
neutron star EOS does not drastically depend on the behavior of E sym at high density. 

A similar analysis for the standard RMF plus an additional isovector-vector nonlinear 


term model of Horowitz and Piekarewicz [1] with Z271* parameter set has been also per¬ 
formed. The various sets of coupling constant combinations are shown in Table m whereas 
their effects on matter properties are shown in Fig [21 Similar conclusion is obtained both for 
proton fraction and M* in the PNM. Significant dependency in the EOS of this model on 
E sym is found. The different trend in E sym and EOS of this model compared to the E-RMF 
one originates from the different form of the isovector-vector non linear terms used in both 
models. 

To investigate which factor dominantly controls the behavior of the NMFP based on 
RMF models, we compare the EOS and NMFP of both models with the one of the standard 
RMF model, i.e., the NL-3 parameter set. The results are shown in Fig|21 In the upper left 
panel, we show the effects from the variation of E sym in the G2* parameter sets (represented 
by the variation of g p and rj p values) on the NMFP. Similarly, for the Z271* and NL-3 
parameter sets (represented by the dots form), the results are shown in the upper right 
panel. The EOS and M* in the PNM of the standard G2, Z271 and NL-3 parameter sets 
are displayed in the lower left and right panels, respectively. Different from NL-3, which has 
an anomalous behavior in its NMFP, it is found that G2* and Z271* parameter sets predict 
the NMFP trends, which do not change with the variation of E sym (in both models, the 
anomalous behavior in the NMFP does not appear in every parameter set even though they 
have different E sym values). This means that the appearance of the anomalous behavior in 
the NMFP seems to be insensitive to the value of the proton fraction. The NL-3 parameter 
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/Wpo 


/Wpo 




e(GeV fm ) 


/W Po 


FIG. 2: Isovector-vector channel adjustment in the Horowitz and Piekarewicz model (Z271) [jj. 
Effects of different combinations of g p and Ay on the symmetry energy of SNM are shown in the 
upper left panel, on pressure and M* of PNM in the lower left and right panels, respectively, and 
on proton fraction predictions in the upper right panel. Shaded region in the upper right panel 
corresponds to the proton fraction threshold for the direct URCA process 


set has PNM with a stiff EOS and a relatively small M* value but, on the contrary, Z271 
and G2 have PNM with a soft EOS and large M* value at high density. This fact gives an 
indication that a soft EOS and a normal behavior of the NMFP are mostly determined by 
the relatively large M* value at high density. On the other hand, finite nuclei calculations 
using the standard RMF model Jjj. []JJ, 120, bid] inform us that acceptable shell structure 
predictions in finite nuclei regions require a small M* value (~ 0.6 M) in the saturation 


density which is fulfilled by G2 


291 and NL-3 


3 it 0 E 


m parameter sets. 


Actually, the above results could be more clearly interpreted by looking at the effects of 
different M* on the NMFP in the same model. The difference in each parameter set is only 
in the values of some adjusted coupling constants, while they should have acceptable SNM 
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/W Po f ) ?J Po 




e(GeV fm 3 ) Pb'Po 


FIG. 3: Effects of the isovector-vector channel adjustment on NMFP predictions for G2, NL-3 
and Z271 models. The upper left panel is for the E-RMF model (G2) and the upper right panel is 
for the Horowitz and Piekarewicz model (Z271). For comparison, the result of NL-3 is given in the 
upper-right panel by the dots form. The pressure as a function of the energy density and M* as a 
function of the nucleon-saturation densities ratio are given in the lower left and lower right panels, 
respectively. 


TABLE V: M* variations in the NL-Z parameter set of the standard RMF model. 


Parameter 

NL-Z 

P-070 

P-080 

9 s/( 4tt) 

0.801 

0.673 

0.575 

gv/( 4vr) 

1.028 

0.817 

0.632 

K-3 

2.084 

2.463 

2.899 

K 4 

-8.804 

-7.595 

12.779 
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e(GeV fm' 3 ) 



Pb^Po 


FIG. 4: Effects of the M* variations on the proton fraction and on the NMFP (upper left and 
upper right panels), and on the EOS (lower left panel). Shaded region in the upper left panel 


corresponds to a threshold of the proton fraction for the direct URCA process 


33]. Variations of 


the M* as a function of the nucleon-saturation densities ratio is shown in the lower right panel. 


and PNM predictions at saturation density (more details about matter predictions of these 
models can be seen in Ref. [ 37 1 ), and for our purpose here, the proton fraction predictions 
should be similar. The parameter sets of Ref. |22| are suitable for this task. The coupling 
constants variations of the models are shown in Table El whereas the results are shown 
in Fig. 01 The effects of the M* variation on the proton fraction are shown in the upper 
left panel, while the effects on the NMFP and EOS of PNM are shown in the upper right 
and lower left panels, respectively. Clearly, if M* becomes too low then the corresponding 
EOS becomes too stiff and the anomalous behavior in the NMFP appears. The P-070 and 
P-080 parameter sets have unacceptable shell structure predictions in some nuclei 131] since 
these parameter sets have a too large M* in the saturation density which, as a consequence, 
leads to a too narrow spin-orbit splitting prediction. Recently, “FSU GOLD” parameter 
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TABLE VI: Effects of the 5 meson on the G2 parameter set of the E-RMF model. Case g p = 0.39. 


Isovector 


Set 


parameter 

I 

II 

III 

9p 

9.483 

12.313 

15.937 

96 

0 

5.026 

7.540 

•, 5 meson in the G2* parameter sets of the E-RIv 

Isovector 


Set 


parameter 

I 

II 

III 

9p 

11.786 

15.304 

18.784 

96 

0 

5.026 

7.540 


set has been introduced by Todd-Rutel and Piekarewicz |JJ| which yields a soft EOS, while 
still accurately reproducing experimental data of binding energies and charge radii of some 
magic nuclei and also centroid energies for breathing mode of 208 Pb and 90 Zr. Unfortunately 
the shell structure prediction of this parameter set is not reported in that paper. Therefore, 
before drawing any further conclusion, we should wait for their full calculation, including 


the predicted shell structure properties of some magic nuc 
In conclusion, these results confirm previous findings 


ei. 


29 


30] about the wide range of 


applications of the E-RMF model. In our view, the reason comes from the fact that the 
E-RMF model has a relatively small M* (~ 0.6 M ) in saturation density (demanding feature 
for finite nuclei) but a relatively large M* at high density (demanding feature for the neutron 
star). Extra nonlinear and tensor terms of this model compared to the standard one seem 
to be the source of this behavior (see Table UJ) . From the possibility that the density depen¬ 
dent E sym can be adjusted, the claim that RMF models predict relatively lower threshold 
densities for the direct URCA process and this fact can be considered as a weak point of 


the models 


38, 


39 


Q, 


now can be re-explored. A precise density dependent E sym at high 


density experimentally determined and/or extracted from the properties of the neutron star 
are needed in this case. 

To study the effects of 5 meson on the E-RMF model, we start with the standard G2 with 
7] p = 0.39 and generate different g p and gs parameters combinations but we keep the same 
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/Wpo 



/Wpo 




e(GeV fm 3 ) Pb'Po 


FIG. 5: Effects of the 5 meson addition in the E-RMF model for g p = 0.39 on the symmetry 
energy and proton fraction (upper left and right panels) and on the EOS and M* in the PNM 
( lower left and right panels). Shaded region in the upper right panel corresponds to the proton 
fraction threshold for the direct URCA process Q- 

E sym = 26.57 MeV at kp = 1.172 fm -1 (note: there would be no change in the conclusion if 
we took E sym = 24.1 MeV at kF = 1-14 fm -1 ). The combinations of the coupling constants 
of the models can be seen in Table eh The matter properties predictions are shown in 
Fig. E It is clearly seen that the presence of the 5 meson results in a higher E sym at high 
density. This fact leads to a higher proton fraction. But, in the region pB=(l-5-2)p 0 , the 
difference between the E sym value of the E-RMF plus a 5 and that without a 5 meson is not 
so significant. The presence of the 5 meson also makes the PNM EOS stiffer and reduces 
the value of the PNM M*. The reduction magnitude depends on the magnitude of gg . A 
similar trend is also found in the case of r] p = 4.49. The combinations of g p and g$ coupling 
constant are shown in Table rvTTl and the corresponding results can be seen in Fig. El 

In Fig.0 it can be seen that the presence of the 6 meson removes the anomalous behavior 
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0.25 

0.2 
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>h 0.15 
0.1 
0.05 

1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



/Wpo 


/Vai 




e(GeV fm 3 ) P^Po 


FIG. 6: Same as in Fig. El but for r\ p = 4.49. 


in the predicted NMFP. The effect appears more pronounce in the case of ij p = 4.49, rather 
than ri p = 0.39. This fact is clearly depicted in Fig. 0 

In Fig. 0 we show the E/A ratio and the pressure, for the SNM as well as the PNM of 
the four RMF parameter sets. NL-3 is a parameter set with good predictions for observables 
of finite nuclei and has a stiff EOS at high density. Z271 is a parameter set that is specially 
constructed for the neutron star and has a soft EOS at high density. G 2 is a parameter set 
with acceptable predictions for observables of finite nuclei and has a relatively soft EOS at 
high density. Clearly, parameter sets with soft EOS are consistent with experimental data 


of Danielewicz et al. 
et al. 


44j and close to the results of the variational calculation by Akmal 


41], BHF calculation with AV14 potential plus 3BF of Baldo et al. j3], and DBHF 


calculation of Li et al. {42j. It seems also from the results of G2* + 5, that the enhancement in 
the isovector channel of G2 shifts the E/A ratio and pressure predictions of that parameter 
set closer to the result of variational calculation from Akmal et al. 141], 
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/W Po 


FIG. 7: Effects of the 6 meson addition in the E-RMF model on the corresponding NMFP 
prediction. Left panel is for tj p = 0.39, while right panel is for rj p = 4.49. 





FIG. 8: Energy per nucleon {E/A) of the SNM (upper left) and PNM (upper right) of G2, Z271, 
NL3 and G2* + 5 parameter sets. The corresponding pressures are given in the lower left and lower 
right_panels. For comparison, we also show the results from variational calculation of Akmal et 
al. 41], DBHF calculation of Li et al. |2]i and BHF with AV14 potential plus 3BF of Baldo et 


al. 


43]. Shaded regions correspond to experimental data from Danielewicz et al 


• Ci- 
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IV. SUMMARY 


In summary, we find that by adjusting the parameters in the isovector-vector sector of 
RMF models we can obtain a low proton fraction at high densities in neutron star. The 
anomalous behavior in the NMFP of RMF models would not appear, if their M* predictions 
at high density were sufficiently large. The presence of the 5 meson in the E-RMF model 
increases the proton fraction at high density but the change in the value of the proton 
fraction is not significant in the region p j g=(1.5-2)po- For the PNM case, the presence of 
the 5 meson has effects that the EOS prediction becomes stiffer and M* becomes smaller at 
high density. The presence of the 6 meson in the E-RMF model also removes the anomalous 
behavior in the NMFP. By adjusting the isovector-vector sector and adding the 6 meson 
in the E-RMF model, the E/A ratio and the predicted pressure of the PNM become much 
closer to the results of Akmal et al. ju]]. 
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